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(N ABSTRACT 
O 

^ We examine a low energy mechanism for the transfer of meteoroids between 

^ two planetary systems embedded in a star cluster using quasi-parabolic orbits of 

^ minimal energy. Using Monte Carlo simulations, we find that the exchange of 

meteoroids could have been significantly more efficient than previously estimated. 

Our study is relevant to astrobiology as it addresses whether life on Earth could 
Oh have been transferred to other planetary systems in the solar system's birth 

^ cluster and whether life on Earth could have been transferred here from beyond 

'p_^ the solar system. In the solar system, the timescale over which solid material 

Q was delivered to the region from where it could be transferred via this mechanism 

^ likely extended to several hundred million years (as indicated by the 3.8-4.0 Ga 

5^ epoch of the Late Heavy Bombardment). This timescale could have overlapped 

with the lifetime of the Solar birth cluster (~ 100-500 Myr). Therefore, we 
^ conclude that lithopanspermia is an open possibility if life had an early start. 

Q\ Adopting parameters from the minimum mass solar nebula, considering a range of 

^ planetesimal size distributions derived from observations of asteroids and Kuiper 

^ Belt Objects and theoretical coagulation models, and taking into account Oort 

^ Cloud formation models, the expected number of bodies with mass > 10 kg that 

^ could have been transferred between the Sun and its nearest cluster neighbor 

^ could be of the order of 10^^-3-10^^, with transfer timescales of 10s Myr. We 

estimate that of the order of 3 ■ 10^ ■ /(km) could potentially be life-bearing, 
^ where / is the depth of the Earth crust in km that was ejected as the result of 

the early bombardment. 
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1. Introduction 

From the collection of thousands of meteorites found on Earth, approximately one hundred have been 
identified as having a Martian origin, and more than 46 kg of rocks have a lunar origin. The study of 
the dynamical evolution of these meteorites agrees well with the cosmic ray exposure time and with their 
frequency of landing on Earth. A handful of meteorites has also been identified on the Moon and Mars 
(e.g. McSween 1976; Schroder ct al. 2008). These findings, together with dynamical simulations (Gladman, 
1997; Dones et al. 1999), indicate that meteorites are exchanged among the terrestrial planets of our solar 
system at a measurable level. Sufficiently large rocks may protect dormant microorganisms from ionizing 
radiation and from the hazards of the impact at landing. Laboratory experiments have confirmed that 
several microorganisms embedded in martian-like rocks could survive under shock pressures similar to those 
suff'ered by martian meteorites upon impact ejection (Stofliler et al. 2007; Horneck et al. 2008). Other 
studies have also shown that microorganisms in a liquid, bacteria and yeast spores can survive impacts with 
shock pressures of the order of GPas (Burchell et al. 2004; Willis ct al. 2006; Hazell et al., 2010; Meyer et al. 
2011). Therefore, it is of interest to consider that the exchange of microorganisms living inside rocks could 
take place among the solar system planets, a phenomenon known as lithopanspermia. Under this scenario, 
life on Earth could potentially spread to other moons and planets within our solar system; conversely life on 
Earth could have an origin elsewhere in our solar system. 

Melosh (2003) investigated quantitatively the probability of transfer taking place between the stars 
in the solar local neighborhood. He found that even though numerical simulations show that up to one- 
third of all the meteorites originating from the terrestrial planets are ejected out of the solar system by 
gravitational encounters with Jupiter and Saturn, the probability of landing on a terrestrial planet of a 

neighboring planetary system is extremely low because of the high relative velocities of the stars and the low 
stellar densities: after the heavy bombardment, he estimates that only one or two rocks originating from the 
surface of one of the terrestrial planets may have been temporarily captured into another planetary system, 
with a 10~^ probability of landing in a terrestrial planet. Therefore, he concluded that lithopanspermia 
among the current solar neighbors is "overwhelmingly unlikely" . 

In a subsequent paper, Adams & Spergel (2005) pointed out that the majority of stars, including the 
Sun, are born in stellar clusters where the probability of transfer would be higher due the larger stellar 
densities and smaller relative velocities compared to those for field stars (including the current solar neigh- 
borhood). The dispersal time of the clusters and timescale for planet formation are comparable: the former 
is approximately T = 2.3(Mciuster/M0)° '^ Myr = (135-535) Myr (for a cluster of 1000-10000 members - 
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Adams 2010), while the latter is of order 100 Myr for terrestrial planets. Therefore, it could be possible that 
rocky material be transferred before the cluster disperses. Adam & Spergel (2005) estimated the probability 
of transfer of meteoroids between planetary systems within a cluster by using Monte Carlo simulations. 
To increase the capture cross-section they assumed that the stars are in binary systems. They found that 
clusters of 30-1000 members could experience O(10^)-O(10^^) capture events among their binary members. 
Adopting typical ejection speeds of ~5 km/s, and the number of rocky ejecta (of mass > 10 kg) per system 
of Nj^ ^ 10^®, they found that the expected number of successful lithopanspermia events per cluster is ~ 
10~^; for lower ejection speeds, ^ 2 km/s, this number is 1-2. These latter estimates are relevant to the 
exchange of biologically active material. Valtonen et al. (2009) have also studied the exchange of solids 
between stars in the solar birth cluster and its enhanced capture probability compared to the exchange of 
solids between field stars. They concluded that approximately 10^*^ bodies with sizes larger than ~ 20 cm 
may have been exchanged between the cluster stars (compared to 10^* between field stars). 

Because there is a significant increase in the number of possible transfer events with decreased ejection 
velocity, it is of interest to study a very low energy mechanism with velocities significantly smaller than 
those considered in Adams & Spergel (2005). This mechanism was described by Belbruno (2004) in the 
mathematical context of a class of nearly parabolic trajectories in the restricted three-body problem. The 
escape velocities of these parabolic-type trajectories are very low (^ 0.1 km/s), substantially smaller than 
the mean relative velocity of stars in the cluster, and the meteoroid escapes the planetary system by slowly 
meandering away. This process of "weak escape" is chaotic in nature. Weak escape is a transitional motion 
between capture and escape. For it to occur, the trajectory of the meteoroid must pass near the largest 
planet in the system. "Weak capture" is the reverse process, when a meteoroid can get captured with low 
velocity by another planetary system. The fact that the escape velocities of the meteoroids we consider 
here are small significantly enhances the probability that a meteoroid can be weakly captured by another 
planetary system, due to the lower approach velocity to the neighboring stars. 

This paper studies the slow chaotic transfer of solid meteoroids between planetary systems within a star 
cluster, focussing on the transfer probabilities which are a critical factor in the assessment of the possibility of 
lithopanspermia. We consider the observed size distributions of asteroids and Kuiper belt objects as well as 
velocity distributions of solar system ejecta from dynamical simulations, to estimate the number of very low 
velocity ejecta, and thereby estimate the number of weak transfer events in the solar birth cluster. Section 
[2] describes the model for minimal energy transfer of meteoroids between two stars in the cluster, where 
the transfer takes places between two "chaotic layers" around each star; these chaotic layers are created by 
the gravitational perturbations from the most massive planet in the planetary system and from the rest of 
the cluster stars. We refer to this transfer mechanism as "weak transfer" . Section [3] describes the location 
of these chaotic layers and uses geometrical considerations to obtain an order-of-magnitude estimate of the 
probability that meteoroids that have weakly escaped a star are weakly captured by its nearest neighbor in 
the cluster. The latter calculation is refined in Section |4] using Monte Carlo simulations. Section [5] calculates 
the number of weak transfer events between the young solar system and the nearest star in the cluster; 
we explore two cases, where the target star is a solar-type and a low-mass star. Finally, in Section [6] we 
summarize our results and discuss the implications for lithopanspermia. 
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2. Minimal energy transfer of solid material within a star cluster 

In this section, we first introduce the concepts of weak capture and weak escape (see Belbruno 2004, 
2007, 2010 for a detailed discussion) and then describe a model of how to construct a minimal energy transfer 
of solid material within a star cluster via weak transfer. 



2.1. General planetary system model 

Consider a general planetary system (5) consisting of a central star (Pi) and a system of iV planets 
{Pi,i = 2, ...,N, with > 3) on co-planar orbits that are approximately circular, where the labeling is not 
reflective of the relative distances from Pi. Assume that the mass of the star (mi) is much larger than the 
masses of any of the planets (mj, i.e. mi » rrii, for i = 2, ...,N) and that the mass of one of the planets 
(P2) is much larger than the masses of all the other planets (i.e. m2 ^ m^, for i = 3, - this condition 
is fulfilled in the case of our solar system where P2 is Jupiter). A meteoroid (Pq) is considered to have a 
negligible mass (mo = 0) with respect to the mass of any of the planets and therefore does not gravitationally 
perturb the circular orbit of P2. Without loss of generality, we view S to consist of Pi and the planet P2, 
moving around Pi in a circular orbit of radius A < 500 AU. The motion of Pq within this system is the 
classical three-dimensional restricted three-body problem, hereafter RP3D. If Pq is constrained to move 
in the plane of motion of P2 and Pi, we have the planar circular restricted three-body problem, hereafter 
RP2D. The differential equations for RP3D are well known (see Belbruno 2004). 

Because we are interested in the possibility of lithopanspermia to and from the solar system, in this 
paper we have assumed that P2 is a Jupiter-mass planet. However, recent observational results by the Kepler 
mission indicate that planetary systems with several Neptune-mass planets maybe more common (Batalha 
et al., in prep.). The weak transfer mechanism described in this paper could also be applied in this latter 
case and in some regards enhanced. We leave the study of these cases for future work. 



2.2. Weak capture and weak escape 

A convenient way to define the capture of Pq with respect to Pi or P2 in RP3D is by using the concept 
of "ballistic capture" . The two-body Kepler energy (Ek) of Pq with respect to one of the bodies P;; (fc = 1, 2) 
is 

1 2 "^fe n\ 

^'= = 2^^-7^' 

where Vk is the velocity magnitude of Pq relative to Pk and is the distance of Pq from P^. The Kcplcr 
energy is a function of the trajectory, Ek{<p{t)), where 4>{t) = (r(t), v(i)) is the solution of RP3D for the 
trajectory of Pq in inertial coordinates, t is time, r = (ri, r2, r^) is the position vector from the center of the 
inertial coordinate system to Pq and v = {vi,V2,V3) is the velocity vector of Pq. Ballistic capture takes place 
when Ek < 0, while ballistic escape occurs at the transition from Ek = and Ek > 0. We are interested when 
Pq changes from hyperbolic motion with respect to P2 to ballistic capture. These trajectories are referred to 
as "weak capture" and Belbruno (2004, 2010) showed they are generally unstable and chaotic. The region 
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around P2 where weak capture occurs in position-velocity space is called a weak stability boundary 




This region around P2 {WSB{P2)) results from the gravitational perturbation of Pi. It can be viewed as a 
location where a particle is tenuously and temporarily captured by P2 : the particle will first move around P2 
for a short time with negative Kepler energy that approaches zero and then increasing above zero yielding to 
hyperbolic escape. Ballistic escape from P2 is referred to as "weak escape" because it occurs near WSB{P2). 

As is evident by numerical integration of Pq around P2, weak capture generally occurs for relatively 
short time spans. For example, in the case of the Earth- Moon system (where /i = ^^^-^^ — 0.012), weak 
capture around the Moon occurs for time spans of days or weeks; in the case of the Sun- Jupiter system, the 
timescale would be months to a few years (Belbruno & Marsden, 1997; Belbruno, 2007). 

In our solar system, the existence of the weak stability boundary and the viability of weak capture was 
demonstrated by the Japanese spacecraft Hiten: using a trajectory designed by Belbruno (1993, 2007), Hiten 
was captured into an orbit around the Moon in 1991 without the use of rockets to slow down. Weak capture 
at the Moon was also achieved in 2004 by the ESA spacecraft SMART 1 (Racca, 2003; Belbruno, 2007). In 
another application, weak escape from the Earth's L4 (or L5) Lagrange point was invoked to suggest a low 
energy transfer to the Earth for the hypothetical Mars-sized impactor that is thought to have triggered the 
"giant impact" origin of the Moon (Belbruno & Gott, 2005). 



In the absence of P2 in the RP2D, a parabolic orbit Q{t) for Pq around Pi is a planar two-body Kep- 
lerian parabolic trajectory. When P2 is considered, there is a "chaotic layer" in Q-space near the parabolic 
trajectory that consists of infinitely many parabolic and near-parabolic trajectories (Xia, 1992; Belbruno, 
2004) . This layer has a positive measure and a two-dimensional transversal slice in four-dimensional position 
and velocity space yields a fractal structure. This layer is obtained when a parabolic trajectory is approx- 
imately at its periapsis with respect to Pi, and where it also has an approximate periapsis passage with 
respect to P2, slightly beyond P2 (i.e. at a radial distance from Pi slightly larger than A). The periapsis 
passage near P2 is done in such a way so that Pq is slightly hyperbolic with respect to P2 (approximately 
in WSB{P2)). This imparts a small gravity assist to Pq. In this case, there are infinitely many possible 
trajectories near the original parabolic trajectory, some which do not move out infinitely far from Pi and 
eventually fall back towards Pi for another possible flyby of P2. Other near-parabolic trajectories will escape 
Pi on hyperbolic trajectories with respect to Pi with a positive escape velocity, a. Because these hyper- 
bolic trajectories lie near parabolic trajectories, a will be small. We refer to these as "low velocity escape 
trajectories" or "minimal energy escape trajectories" . Within the chaotic layer mentioned above exist in- 
finitely many parabolic trajectories that hyperbolically escape Pi with small escape velocity cr. The distance 



P = Pcsc('7^l) from Pi at which Pq achieves this escape velocity can be obtained by the formula a ~ y ? 
where G is the gravitational constant. The fact that these escape trajectories are nearly parabolic and have 
very small escape velocities implies that their Kepler energy with respect to Pi is nearly zero, for sufficiently 



SA detailed mathematical explanation of the WSB can be found in Belbruno et al (2010) and Garcia & 
Gomez (2007). They showed that the weak capture boundary is a complicated region that has a fractional 
dimension and that is equivalent to a Cantor set. The WSB region can be viewed as a limit set of the stable 
manifolds to the Lyapunov orbits associated to the coUinear Lagrange points Pi, L2 of P2. 



2.3. Parabolic motion and low velocity escape from S 
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large R. 

In summary, in the circular restricted three-body problem (Pqj Pii P2), there is a chaotic layer that 
replaces the regular parabolic trajectories of the two-body problem {Pq, Pi); weak transfer can take place in 
this chaotic layer because this set of trajectories have low escape velocities. 



2.4. Model for the minimal energy transfer of solid material within a cluster 

To achieve a minimal energy transfer of the meteoroid Pq from S to another system S* , consider that 
the system S is embedded in a star cluster; this will introduce additional gravitational perturbations on Pq 
which will interact with the gravitational field of Pi and form a weak stability boundary around Pi where the 
motion of Pq will be chaotic. Because the transfer mechanism requires low relative velocities, we consider an 
open star cluster with a low velocity dispersion, with relative stellar velocities U ~ 1 km/s (for comparison, 
older globular clusters have stellar velocity dispersions of several tens of km/s - Binney & Tremaine 1988, 
Meylan & Heggie 1997). After the cluster disperses, the relative distances and relative velocities between 
a star and its closet neighbors increase significantly. For example, the Sun's current closest neighbor a- 
Centauri (not necessarily a cluster sibling) is now 2.6 • 10^ AU away (4.28 pc) moving with a relative velocity 
of 6 km/s. The latter is significantly higher than the ^ 1 km/s required for the weak transfer mechanism. 
In this paper we consider the transfer of solid material between stars in a cluster before the cluster starts to 
disperse. 

Imagine that Pq makes a minimal energy escape from Pi in the Pi , P2 orbital plane by making a slightly 
hyperbolic flyby of P2 at WSB{P2) (where the latter is the weak stability boundary around P2 created by 
the gravitational perturbation of Pi); Pq moves away from Pi with a escape velocity a near zero. Pq then 
gets to WSB{Pi), the weak stability boundary of Pi caused by the gravitational perturbation of the other 
— 1 stars in the cluster (where C^n-i represents this set of N-1 stars, {P2, P3, Pat}); this boundary is 
located at a distance Rcsd'm-i) from Pi, at which the combined gravitational force of the stars in the cluster 
is comparable to the gravitational force of Pi . As described previously, the motion in this region is chaotic 
and lies in the transition between capture and escape from P^ Although WSB{Pi) is a complicated fractal 
set, non-spherical in shape, the sphere of radius RcsciiT^i) around the central star lies approximately in the 
part of WSB{Pi) where the motion is least stable and the trajectories are slightly hyperbolic with respect to 
Pi (Belbruno et al. 2010; Garcia & Gomez, 2007) . We refer to trajectories that escape Pi from ^^^^(Pi) 
as "weak escape" trajectories (note that low-velocity escape does not imply weak escape unless it occurs at 
a weak stability boundary). The structure of the set of near-parabolic trajectories that escape Pi by flyby 



'^The sensitivity of the motion of Pq at the distance R can be deduced from an analogous four-body 
problem described in Belbruno (2004) and Marsden & Ross (2006) for a transfer to the Moon used by the 
spacecraft Hiten. In this case, the four bodies are the Earth (Pi), Moon (P2), spacecraft (Pq), and the 
Sun, analogous to C5n-i. The spacecraft leaves the Earth, flys by the Moon, where the flyby is weakly 
hyperbolic, and then travels out to roughly 1.5T0^ km where the gravitational force of the Sun acting on 
the spacecraft approximately balances that of the Earth. At this location, the motion of the spacecraft is 
highly sensitive to small differences in velocity and lies between capture and escape from the Earth, i.e. lies 
at the weak stability boundary between the Earth and Sun. The motion of Pq in this region is chaotic and 
the path of Pq can be altered using a very small amount of fuel. 
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of P2 yield infinitely many such trajectories. In particular, this implies that for any point p on the circle of 
radius i?esc(wi) around Pi, there is a nearly parabolic trajectory that will pass by p with velocity a, weakly 
escaping Pi. This is illustrated in Figure 1 (see also Figure 1 in Moro-Marti'n & Malhotra, 2005). 

We now compute the trajectory of Pq from the distance WSB{Pi) until it encounters another planetary 
system S* centered on star of mass m\. This is done by using a different set of differential equations; 
during this motion the trajectory of Pq is relatively undisturbed. The meteoroid is then weakly captured 
at WSB{Pl) at the distance R* — Rcapiml) from P^*. Figure 1 gives a schematic representation of this 
process. After Pq is weakly captured at the distance Pcap('7ii) from Pj*, say at time t = T, it moves towards 
Pi for t > T. Analogous to system S, we assume that system S* has a dominant planet P2* orbiting Pj* in a 
circular orbit at a radial distance A*. As t increases from P, Pq will move to a periapsis distance r* from P^*, 
where r* can range from (collision with Pj*) to Rcupiml). However, unlike the escape of Pq from Pi, where 
Pq had a periapsis near the location of P2 and in the same orbital plane of P2 , the approach of Pq into the 
S* system is three-dimensional and need not be restricted to lie in the P| orbital plane, nor have a periapsis 
of Pj* near P|. The resulting motion of Pq for t > T will, in general, be complicated and a priori not known. 
There is no way to predict if Pq will pass near or collide with Pj* without doing numerical simulations that 
are outside the scope of this paper. For this reason, we conclude our search at t = T, when Pq is weakly 
captured at Pi , and calculate the probability of this weak capture to take place. An order-of-magnitude 
estimate of the weak capture probability is derived in Section [3] based on geometrical considerations, while 
in Section |4] we calculate the probability numerically using Monte Carlo simulations. 



3. Order of magnitude estimate of the weak capture probability 

3.1. Approximate location of the weak stability boundary 

We calculate Pcsc("^i) and Pcap(™i) as a function of the stellar mass. For weak escape to take place, 
the velocity a of the meteoroid at the distance Pesc(wi) from the star Pi must be sufficiently small. Because 
we are considering slow transfer within an open cluster with a characteristic dispersion velocity J7 « 1 km/s, 
we require that a is significantly smaller than U , i.e. of the order of 0.1 km/s. This is much smaller than the 
nominal values of several km/s used by the Monte Carlo methods in Melosh (2003) and Adams & Spergel 
(2005). 

To place the above choice for a in context, we study the velocity distribution of weakly escaping test 
particles from the solar system using a three-body problem between the Sun (Pi), Jupiter (P2) and a massless 
particle (Pq). To be consistent with our framework, we model this using RP2D, where the test particle moves 
in the same plane of motion as Jupiter, assumed to be in a circular orbit at 5 AU. The trajectory of the test 
particle is numerically integrated by a standard Runge-Kutta scheme of order six and numerical accuracy of 
10~* in the scaled coordinates. The initial conditions of the test particle is an elliptic trajectory very close to 
parabolic with periapsis distance rp=5 AU and apoapsis distance ra=40,000 AU. (Note that such orbits are 
not dissimilar to those of known long period comets in the solar system.) For each numerical integration, we 
assume that Jupiter is at a random point in its orbit when the test particle starts from apoapsis at 40,000 
AU and falls towards Pi . We record the velocity of the test particle at the time it escapes with respect to the 
Sun (i.e. when the Kepler energy with respect to the Sun is positive) after performing a sufficient number of 
Jupiter fly-bys. Figure 2 shows the distribution of Voo'- out of 670 cases, 58% have Voo < 0.1 km/s and 79% 
have Woo < 0.3 km/s. Based on these results, we will assume the velocity a of the meteoroid at the distance 
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Rcsc{fni) from the star to be in the range 0.1-0.3 km/s. 

For a given a, the location of the part of the weak stability boundary that gives the most unstable 
motion is approximately given by the sphere of radius Rcsdnii) = 2Gmi/a^ , where mi is the mass of the 
star and a is in the range 0.1-0.3 km/s (see Figure 3). Beyond this boundary, when Pq weakly escapes Pi, 
we assume that the meteoroid will move at a nearly constant velocity a with respect to the star. 

Slow transfer to a neighboring planetary system enables the meteoroid to arrive at the distance -Rcap('7ti) 
with a relative velocity with respect to the target star (Pj* ) that is similar or smaller than its parabolic escape 
velocity at that distance, a* = \J~^^^^^^ where m\ is the mass of the target star; if its velocity is higher 
than cr*, it will not be captured and will fly by. Since the relative velocity between stars in the cluster is 
U ~ \ km/s, the meteoroid that weakly escaped from star Pi moves toward the target star P^ with velocity 
U ± a (see Figure 1). Because a is small relative to U and can be neglected, the relative velocity of the 
meteoroid with respect to the target star is « U . Therefore, weak capture can occur at the distance at 
which cr* w [/ w 1 km/s, i.e. i?cap(TO^) = Fi gure 3 shows Rcsd'mi) and RcapimX) as a function 

of the stellar mass. The horizontal lines indicate the range of cluster sizes (dashed-dotted lines) and mean 
interstellar distances (dotted lines) for clusters consisting of N = 100, 1000 and 4300 members respectively; 
clusters with this range of sizes are the birthplaces of a large fraction of stars in the Galaxy (Lada & Lada 
2003). The radius of the cluster depends on the number of stars and is given by i?ciustcr — lpc^N/300 
(Adams 2010). For N=100, 1000 and 4300, we get that i?ciuster is approximately 2.1 • 10^ AU, 6.5 • 10'"^ AU 
and 1.3 • 10^ AU respectively. We can estimate the average interstellar distance within a cluster, D = ■nT^I'^ , 
where n = ^^^^ — is the average number density of stars in the cluster. D is approximately 7T0^ AU, 10^ 
AU and 1.2-10'^° AU, for a cluster with N=100, 1000 and 4300 members respectively. The choice of N =4300 
is explained in Section |4j2j 

It is of interest to compare the results in Figure 3 to those in Melosh (2003). The latter considers the 
exchange of meteoroids between field stars using hyperbolic trajectories; it estimates a cross-section of 0.025 
AU^ for meteoroid capture by a planetary system with a Jupiter-mass planet located at 5 AU. Under the 
scenario considered in the present work, the order-of-magnitude estimate for the weak capture cross-section 
would be TT • i?cap; for ^ solar-mass star, i?cap ^ AU (Figure 3) so that the weak capture cross-section 

is ^ 10^ AU^ (i.e. many order of magnitudes larger than Melosh 2003 estimate). This indicates that weak 
escape and capture within an open cluster can enhance drastically the probability of transfer. A recent 
study in which enhanced capture probability within a stellar cluster has been invoked is that of Levison 
et al. (2010). They modeled the exchange of comets within a stellar cluster assuming that each star is 
surrounded by a large disk of comets with g = 30 AU and a = 1000-5000 AU, and with a population similar 
to that of the Sun; they conclude that up to 90% of comets in the Oort Cloud may be have an extrasolar 
origin. 



3.2. Constraints on Stellar Masses for Weak Transfer 

From Figure 3, we can set constraints on the stellar mass mi that could allow weak escape from Pi to 
take place. The idea is simple: if for a given a (within the range 0.1-0.3 km/s, see Section 3.1 1, we have that 
i?csc('77.i) < -D, i.e. the weak escape boundary is located within the distance of the next neighboring star 
Pi , then weak transfer is possible because at the time the meteoroid passes near the star P^ , its velocity is 
similar to the mean stellar velocity dispersion (/7 1 km/s) and there is a significant probability of capture 
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(which we quantify in Section 3.3 1. Conversely, weak transfer by the process schematically represented in 
Figure 1 is much less likely to take place if i?osc("T.i) > D because at the time the meteoroid passes near 
a neighboring star, its velocity is too high and, as a consequence, it will likely fly by. Figure 3 shows that 
for a — 0.1 km/s, the condition i?esc('7ii) < D for weak escape is satisfied for mi < O.4M0, mi < O.57M0 
and mi < O.68Af0, for clusters with N = 100, 1000 and 4300 members respectively. If we consider the 
higher but still acceptable value a = 0.3 km/s, Figure 3 shows that the stellar mass limits for weak escape 
are mi < 3.5M0, toi < S.IM© and mi < 6.IM0, for clusters with N =100, 1000 and 4300 members 
respectively. A meteoroid escaping nearly parabolically from any star with a mass larger than these limits 
will achieve a velocity of cr < 0.1-0.3 km/s only at a distance larger than the mean interstellar distance in 
the cluster; weak transfer is not likely to take place under such conditions. 

Of particular interest is the case of the solar system as the source of the meteoroids. It has been 
estimated that the Sun's birth cluster consisted of N=4300±2800 members (Adams 2010). For a lAf© star 
in such a cluster, we find that the parabolic escape velocity at the mean interstellar distance of D = 1.2-10^ 
AU is cr = \/2Gmi/D — 0.12 km/s. Because this values of a lies within the range of values of interest 
for weak escape, we conclude that the meteoroids originating in the early solar system could have met the 
conditions for weak escape in the Sun's birth cluster. 

We have established the range of stellar masses that could in principle allow weak escape to take place. 
We now estimate the probability of weak capture by a neighboring planetary system. Whether or not the 
transfer of meteoroids from one star to any star of a given mass mi takes place will depend on the relative 
capture cross-section of the target star and on the number of potential targets. We discuss these factors 
below. 



3.3. Relative capture cross-section of the target star 

In the case of weak transfer, assuming that the meteoroids escape isotropically from Pi, the relative 
capture cross-section of the target star would be given by Cs — G f{Rcap{ra\) / , where D is the distance 
between the two stars, i?cap is the weak capture boundary (illustrated in Figure 1 and estimated in Figure 
3) and the factor G f accounts for gravitational focussing. The assumption of isotropy is not justified in the 
framework of the restricted three-body problem, but it is expected in the presence of the cluster stars and of 
the galactic potential that can isotropize the escaping meteoroids' orbits, in analogy with the isotropization 
of the Oort Cloud of comets (Dones et al. 2004). 

The gravitational focussing factor is given by 

Gf = l + C—)\ (2) 

Voo 

where is the velocity at infinity and Vesc is the escape velocity at the distance i?cap("^i) from Pi*. In 
the case of weak capture, the term G/ increases the cross-section due to enhanced gravitational focussing. 
For example, as we saw previously in the case of the Sun and Jupiter, Vao ~ 0.1 — 0.3 km/s, while at the 
distance Pcap = 40000 AU we have that Ucsc ~ 0.2 km/s. In this case we get that G/ « 2, doubling the 
capture cross-section. This situation occurs in our capture methodology. In the formulation for determining 
Rciip{'m-i)i the meteoroid has an approximate relative approach velocity to the target star Pi of roughly 
U = I km/s, which represents Wqo. However, Pcap is determined so that this same value of velocity is taken 
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for Vcsc from the target star. That is, we are assuming v^o ~ v^sc ~ 1 km/s. This imphes also that G/ — 2. 
This is a conservative estimate and does not make use of the nature of weak capture dynamics. In fact, the 
value of G/ can be substantially increased, as we demonstrate below, if at a given value of Rca.p{ml), Voo 
is smaller than the approximate value of U — 1 km/s, while i;esc ~ C/; in this case, — )• which implies 
G/ cx). 



3.4. Number of potential targets 

We denote PimfCw*) as the probability of finding a star of mass ml in the star cluster. This probability 
is referred to as the initial mass function (IMF) and it can be inferred from observations of stellar clusters. It 
is found that a wide range of clusters, varying from large clusters like the Trapezium to smaller clusters like 
Taurus, as well as older field stars, show very similar distributions of stellar masses down to the hydrogen 
burning limit at ~ 0.1 Mq (Lada & Lada 2003 and references therein). The IMF can be characterized 
by the broken power-law, ^(M) = ^iM'^ ^ for 0.6-100 Mq, and ^(M) = C2M-^-i 0.1-0.6 Mq, where 
f(M)dM is the number of stars with mass (M, M+dM). There is a steep decline into the substellar brown 
dwarf regime and a possible second peak but we will ignore objects below the hydrogen burning limit. The 
resulting average stellar mass is ~0.88 Mq. To calculate Pimf(?71i) (square symbols in Figure 4), we use a 
logarithmic binning of masses with d(logM)=0.1 and normalize the distribution to unity, which gives ^i=0.19 
and ^2=0.34. 

An upper limit to the probability that a meteoroid escaping from a star of mass mi (within the range 



described in Section 3.2 ) will get captured by a neighboring star of mass ml is approximately given by the 
relative capture cross-section of its weak stability boundary, G/(i?cap(wi)/-D„i»_„jj)^, where i?cap is the 
radius of the weak stability boundary for capture of the target star and Dmi-mi is the average distance 
between two stars of masses m\ and mi respectively. If we assume that this distance is the average interstellar 



distance D estimated in Section 3.1 we get that for a solar- type star, an upper limit to the weak capture 
probability is approximately 1.3•10-^ 6.5-10-'' and 4.5-10-'', for a cluster of N = 100, 1000 and 4300 
members respectively (see dashed line of Figure 4). However, the distance to a star of a given stellar mass 
will not be D necessarily, but it will depend on the distribution of stellar masses in the cluster (the IMF). 
The simplest case is when both stars have equal masses mi — mi. In this case, the average interstellar 
distance would be ~ (l/'^mi)^^'^! where rimj is the average number density of stars with mass mi, 

37V 

= A ga"'^ — ' and Nm. is the total number of stars in the cluster with mass mi, Nm, = N ■ Pimf(wi). 

-"-cluster 

Because — n • ^imf(™i), we get that ^ D ■ {PiMF{mi))'^^^^ . This means that the transfer 

probability between two stars of equal mass is given by G f{Rca.p{ml)/ D)^ ■ (Pimf("^i))^/'^, where D is the 
average distance between any two stars in the cluster (regardless of their mass). As we have discussed above, 
the value of the focussing factor Gf can be very large. A conservative estimate of the capture probability can 
be done by assuming G/ = 2. Figure 4 shows that the capture probability between two planetary systems 
with solar-type central stars {ml = nii = 1 Mq) are LO-IO-"*, 8.1-10-^ and 5.6-10-^, for a cluster of N = 
100, 1000 and 4300 members respectively. 
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4. Numerical estimate of the weak capture probability based on Monte Carlo 

simulations 

The order-of-magnitude estimate of the weak transfer probabihty discussed in Section[3]has two impor- 
tant caveats: (a) it assumes that the capture takes place in the orbital plane of Pj* (the primary planet in the 
capturing system), i.e. it does not consider the 3-D problem; and (b) the adoption of a focussing factor G/ 
= 2 is conservative because it does not make use of the nature of nonlinear weak capture dynamics, where 
the gravitational forces of Pi, , and the stars of the cluster are all acting on Pg- Gf can increase very 
significantly if at a given value of i?cap(TOi), Voo is much smaller than the approximate value of C/ = 1 km/s, 
and can become 0, while Uosc ~ U (for an example, see Appendix A). In this section, we refine the estimate 
of the weak transfer probability using Monte Carlo simulations that address the caveats mentioned above 
by considering the more general and realistic model that calculates the motion of Pq, Pi, Pj* by a general 
three-dimensional Newtonian three-body problem plus the effective gravitational perturbation of the cluster 
stars. We sample millions of trajectories for the meteoroid Pq weakly escaping Pi testing whether they are 
weakly captured near Pj* . 



4.1. Modeling procedure 

To model the motion of Pq from the distance Posc('™i) from Pi to S* and with initial velocity cr, we 
consider a general Newtonian three-dimensional three-body problem that is perturbed by the gravitational 
force of the stellar cluster. This model gives the motion of Po,Pi,Pi in an inertial coordinate system 
(ri,r2,r3). We assume the gravitational perturbing force of the cluster is obtained from a spherically 
symmetric Hernquist potential 

Uc{r) = — — , (3) 
r ~\- a 

(Hernquist 1990), where M is the total cluster mass, = + r2 + r^, and a is the cluster scale length, 
which is approximately the radius of the cluster, Pciustcr- The center of mass of the cluster is at the origin, 
r = (0,0,0). 

The system of differential equations modeling the motion of Pg , Pi , Pj* , of mass tuq , mi , TOi respectively, 
is given by 

r/c =Vfc, Vfe = mj^ — , (4) 

where = (rfei, rfc2, J^fcs) and = (wfei, 1^/02, Wfcs) = (»^/ci, rfc2, J^fcs) are the position and velocity vectors, 
respectively, of the fcth particle, Pk (with k = 0, 1, 2; P2 represents Pi and m2 represents ml). The 
gravitational potential is V = U + mkUc{rk), where is the distance of Pk to the origin, = \rk\, and 

c/ = y^^. (5) 

U — C/(ri, r2, r^) is a function of 9 variables r^j, j — 1, 2, 3, and 



dV 
dvk 



( 



dV dV dV 



) 



(6) 



where r^y^. is the distance between Pj and Pfc, r^/j, = |rj — Vk\. 

rriQ is approximately relative to mi, m^. 

In the Monte Carlo model, the trajectory of motion for Pq is determined from Q by providing its 
initial position at the distance R — Rcsci'm-i) from Pi at t = 0. The initial velocity of Pq with respect to 
Pi at the distance R is chosen from a distribution of weak escape velocities (see Figure 2). Pq then weakly 
escapes Pi. The initial positions and velocities of Pi and Pi are also given at t = 0. Their initial separation 
distance depends on the cluster properties and is a function of the number of stars in the cluster N. After 
the initial positions and velocities of Po,Pi,Pi are given, the trajectory of Pq is propagated for t > 0. We 
search for the condition of the first weak capture of Pq with respect to Pi, given by a negative value of 
the Kepler energy of Pq with respect to Pj*, < 0. This implies that for Pq, Voo , and therefore 
Gf — ?> oo. This occurs at the weak stability boundary of Pi, WSB{Pi), at a distance R* — Rcapiml) 
from Pj*. This stability boundary is formed around Pi due to the resultant gravitational perturbation of 
the iV — 1 remaining stars of the cluster. We are interested in cases where the time of propagation T is 
smaller than the cluster dispersal timescale, which is a function of the number of stars in the cluster, N. 
Under the conditions described above and using a Monte Carlo approach, we sample millions of trajectories 
of Pq weakly escaping Pi using Runge-Kutta-Nystrom 12th/10th order, variable step, symplectic integrator 
(Dormand et al. 1987). The Monte Carlo method calculates the number of these particles to be weakly 
captured by Pi, i.e. the probability of weak transfer. 



We adopt those cluster properties thought to be representative of the Sun's birth cluster. These prop- 
erties are inferred from a wide range of physical considerations, including the effect of supernova explosions 
on the enrichment of short-lived radioactive isotopes in the solar nebula, protoplanetary disk truncation 
due to photoevaporation from nearby hot stars, and the dynamical disruption of planetary orbits due to 
close encounters with cluster stars. The observed solar system properties that depend on the Solar birth 
environment (such as the evidence of short-lived radio-nuclides in meteorites and the dynamical properties 
of the outer solar system planets and Kuiper belt) led Adams (2010) to conclude that the Sun was born in a 
moderately large cluster with N = 1000-10000 and (N) ~ 4300 ± 2800 members, similar to the Trapezium 
cluster in Orion. For reference, Adams (2010) notes that approximately 50% of stars are born in systems 
with N > 1000, but for ~ 80% of these stars the clusters dissolve quickly after 10 Myr; only ~ 10% of the 
total number of stars would be born in open clusters with lifetimes of the order of 100-500 Myr. 

Using the initial mass function in Section |3.4( we estimate that the total cluster mass is Mdustcr = 
(3784 ± 2500)Mq. Following Adams (2010), the properties of such a cluster would be as follows. The scale 
length of the cluster is approximately its size, given by a = Pciustor = lpc-\/N/300 — 3.78 ± 1.5 pc. With the 
stellar number density, n = 3N/{4t:R^i^^^^j.) = 19 pc~^, the average interstellar distance is D — n^^/^ = 0.375 
pc. The typical expected distance of the Sun to the center of mass of the cluster would be dcm = 2Pciustcr/3 
— 2.52±1 pc (using a radial profile for stellar density consistent with a Hernquist potential). The cluster 
lifetime is approximately T = 2.3(Mciustor/M0)°-^ Myr = (135-535) Myr, for N = 1000-10000. 



4.2. Cluster properties 
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4.3. Modeling assumptions 

For the Monte Carlo simulations we adopt the average values of the solar system birth cluster: N ~ 
4300, D — 0.375 pc, a — 3.78 pc, dcm — 2.52 pc, and a cluster lifetime T — 322.5 Myr. We also make the 
following assumptions: 

• For the duration of the simulation, the cluster size and therefore the average distance between the 
stars, D, are kept constant (i.e. we making the simplifying assumption that the cluster disperses 
instantly at the end of the simulation, at time equal to the average cluster age). 

• At t = 0, the initial separation of Pi,P* is D (coinciding with the average separation of the cluster 
stars); the center of mass of Pi, Pi, placed on the Qi,Q2 plane for convenience, is assumed to be a 
distance c?cm from the inertial frame origin, (5 = 0= (0, 0, 0). 

• Let vi and be the initial velocity vectors of Pi and Pi respectively (relative to Q = 0). At t = 0, 
these two vectors lie on the plane a^Q and we assume |vi — v^fj = 1 km/s, |vi| < 2 km/s and |v^| < 
2 km/s, with the initial angle 9 between vi,v^ varying over [0, 27r] on a^^. 

• At t — Q, Pq escapes from Pi in three-dimensions with the tangential velocity vector Vq/i (relative 
to Pi); the magnitude of Vq/^ is chosen from the weak escape Uqo distribution (D^^) in Figure 2; 
note that this figure assumed Pi has a mass mi = IM© and a Jupiter size planet orbiting at 5 AU, 
but we also assume this distribution for the case toi = O.5Af0 considered below. At i = 0, when Pq 
escapes Pi, the distance between the two is tq/i — |ro/i| G [0.02,0.2] pc. Let (j)o,do be the spherical 
angles that specify the position of the vector rg/i of Pq with respect Pi (note that 9o is distinct from 
9 and should not be confused with it). We assume (po G [0, 27r] and G [0,7r], where 0o is uniformly 
distributed in [0, 27r], and 9o is sinusoidally distributed over [0,7r]. For i > 0, as Pq escapes Pi in 
three-dimensions, it can approach Pi from any direction. 

• We search for weak capture with respect to Pi in the time interval t £ [0, TmaxlMyr, where T^ax is 
the estimated age of the cluster Tmax = T = 322.5 Myr. The condition for capture is a negative value 
of the two-body Kepler energy of Pq with respect to Pi , < 0. This can occur at any distance from 
Pi , or even at the initial time t = when Pq escapes Pi (the latter can occur because Pq escapes Pi 
between 0.02 pc and 0.2 pc, but Pi is only initially at a distance D from Pi). When weak capture 
occurs, we have obtained a weak transfer of Pq from Pi to Pi*. 

In summary, the Monte Carlo simulation is done by randomly choosing values of the paramaters 
^'o/i,''o/i,^o>o,^', |vi| and |vj| at < = that satisfy: Vq/i G D^^, tq/i e [0.02,0.2] pc, 0o S [0,27r], 
6*0 e [0,7r], 9 € [0,2tt], |vi| < 2 km/s, |vj| < 2 km/s and |vi - v^| = 1 km/s respectively. The Monte Carlo 
algorithm searches for the condition < ioi t G [0, Tmax] Myr. 



4.4. Results from the Monte Carlo simulations 

The Monte Carlo simulations study the capture probability within the cluster described in Section |4.2| for 
the following three different combinations of stellar masses: mi = 1 Mq, ml = 1 Mq (Case 1); mi — 1 
Mq, ml = 0.5 Mq (Case 2); and mi = 0.5 Mq, to^ = 1 Mq (Case 3). These simulations explore the weak 
transfer of material between two solar-type stars and between a solar-type star and a star half its mass, the 
interest of the latter being that low-mass stars are the most abundant in the galaxy, with ^ 64% and 72% 
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having mass < 0.4 M0 and 0.6 Mq respectively. For each case, the Monte Carlo simulation ran 5 million 
trajectories, a sufficient number so that the randomization of the initial values produces distributions of the 
parameters that span their respective ranges. The resulting distributions are discussed in the Appendix B. 

For Case 1 (weak transfer of meteoroids between two solar- type stars), the Monte Carlo simulations 
give a revised estimate of the weak capture probability of 1.5-10~^ (see Table 1). Note that the Monte 
Carlo simulations assume that Pj* with ml = 1 M0 is at the average interstellar distance D = rT^^^ (i.e. 
without regard to the distribution of stellar masses in the cluster). The Monte Carlo derived weak capture 



probability is approximately a factor of 3 larger than the order-of magnitude-estimate in Section 3.4 which 
was approximated by the relative capture cross-section of weak stability boundary of the target star, given 
by G f{Rcap{iTLl)/ D)'^ = 4.5-10^'', for = 1, N = 4300 and using a conservative gravitational focussing 
factor of G/ — 2. This increase in the weak capture probability estimated by the Monte Carlo simulations 
is likely due to the larger focussing factors that result from the nature of weak capture. Note that both 
estimates assume that the meteoroids are initially on weakly escaping trajectories from Pi. For Cases 2 and 
3, the Monte Carlo simulations give a weak capture probability of 0.5-10^'^ and 1.2-10^'^ respectively (see 
Table 1). 



5. Estimate of the number of weak transfer events for the young solar system 

With the improved estimate of the probability of weak capture obtained from the Monte Carlo simula- 
tions, we now calculate the total number of meteoroids that could get transferred between two neighboring 
solar- type stars harboring planetary systems by multiplying the capture probability in Table 1 (for Case 1) 
with the estimated number of meteoroids (Nj^) that a planetary system may eject before the cluster dis- 
perses. Because the scenario considered here assumes that the escaping meteoroids are on weakly escaping 
trajectories, Nr is not the total number of escaping meteoroids but only the subset on weakly escaping 
trajectories. The main uncertainties in assessing whether weak transfer is a viable method for meteoroid 
transfer lie in estimating this number. While this number is highly uncertain, we proceed with estimating 
Nj^ for our solar system, as it is the only planetary system for which observations and dynamical models 
enable us to make an educated estimate. 

We estimate Nr from Oort cloud formation models because Oort cloud comets are weakly bound to the 
solar system and are therefore representative of a population of meteoroids that may have been delivered on 
its weak stability boundary, subject to weak escape. We use the model by Brasser et al. (2011) in which ~ 1% 
of the planetesimals in the Jupiter-Saturn region (4-12 AU) became part of the primordial Oort Cloud in the 
early history of the solar system. Under this scenario, the forming giant planets scattered the planetesimals 
in this region out to large distances where they were subject to the slowly changing gravitational potential of 
the cluster; the latter caused the perihelion distances of the scattered planetesimals to be lifted to distances 
3> 10 AU, where the planetesimals were no longer subject to further scattering events but were also safe 
from complete ejection, and thus remained weakly bound to the solar system. To calculate how many 
planetesimals formed in the 4-12 AU region, we adopt the minimum mass solar nebula (MMSN), with a dust 
surface density E = Eoia/AOAU)-^/'^ , where Eq = Eq^ = 0.1 g/cm^ (Weidenschilhng 1977; Hayashi 1981). 
Integrating between 4 and 12 AU, we find a total mass in solids of 10^^ g. With the above estimate for the 
total mass in solids, we calculate the number of planetesimals by adopting a planetesimal size distribution 
function representative of the early solar system. This size distribution is uncertain, but can be constrained 
roughly from observations and coagulation models. Based on these studies, we adopt three power-law size 
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distribution functions for our calculation of Nr, with dN/dD cx D for _D > Dq and dN/dD cx D ''^ if 
D<Do, of different power-law indices at the small and large sizes: 

• Case A: qi — 4.3, 52 ~ 3.5, Dq = 100 km, Dmax = 2000 km Pluto's size), Dmin — 1 /^ni dust 
blow-out size). This size distribution has the power law index of a collisional cascade at the small 
size end, and that of the hot Kuiper Belt at the large size end (from Bernstein et al. 2004), with a 
break diameter consistent with that of the present-day Kuiper belt. The approximate scenario of this 
case is as follows: in the Jupiter-Saturn zone, the accretion of large planetesimals proceeded to make 
the large-size bodies similar to those observed in the present-day Kuiper belt, whereas at the small 
size end, the dynamical stirring by the large bodies produced a classical collisional cascade. For the 
latter, the reason why we do not adopt the observed present-day power law index of the hot Kuiper 
belt at sizes < 50 km ((72 = 2.8, Bernstein et al. 2004) is because these smaller bodies are the result 
of an advanced erosional process that changed their size distribution on gigayear timescales (Pan & 
Sari 2005). 

• Case B: qi = 3.3, 52 = 3.5, Dq = 2 km, D^ax = 2000 km (^ Pluto's size), Dmin = 1 (~ dust 
blow-out size). This size distribution is derived from theoretical coagulation models (see review in 
Kenyon et al. 2008). 

• Case C: qi — 2.7, q2 — 3.5, i^o = 2 km, D,„ax = 2000 km (^ Pluto's size), D,„in = 1 /ini (^ dust 
blow-out size). This size distribution is also based on theoretical coagulation models (see review in 
Kenyon et al. 2008). 

• Case D: qi — 4.3, q2 = 1.1, Dq = 100 km, Dmax — 2000 km (^ Pluto's size), 'Drnin = 1 ("^ 
dust blow-out size). This size distribution represents perhaps the "worst case" scenario in which the 
depletion of the small bodies took place very early on, before the objects were transferred into the 
Oort Cloud. The parameters here are within the range that Bernstein et al. (2004) find for the Kuiper 
belt. This size distribution is also similar to that discussed in models for the primordial asteroid belt 
{qi = 4.5, (72 = 1-2, Do ~ 100 km, Bottke et al. 2005). 

Of particular interest are the meteoroids >10 kg, that may be large enough to shield potential biological 
material from the hazards of radiation in deep space and from the impact on the surface of a terrestrial planet 
(Horneck 1993; Nicholson et al. 2000; Benardini et al. 2003; Melosh 2003). 

Given a total mass of 10^^ g of solids in the 4-12 AU zone in the solar system at the time of Jupiter 
and Saturn formation, and adopting the above size distributions of the solid bodies, we can calculate the 
number of planetesimals with masses > 10 kg (diameter D > 26 cm if p = 1 g/cm^^). Then, assuming that 
1% of these planetesimals were delivered to the weak stability boundary of the solar system, we obtain an 
estimate for Nr. Table 2 shows these results for the four size distributions considered. We find that Nr is 
in the range 8 • 101^-2 • 10^^ for Cases A, B and C, but is only 5 • 10^ for Case D. 

Finally, multiplying Nr by the weak capture probability determined in Section [4] (1.5-10^'^ for the 
transfer between two solar- type stars), we estimate the number of weak transfer events, A'wte, that could 
have occurred between the early solar system and a neighboring star in the cluster, assuming it is also a 
solar- type star and harbors a planetary system. The results are listed in Table 2; for Cases A, B and C they 
are of the order of 10^^-3 • 10^^, and for Case D it is of the order of 10**. If the target system is a low-mass 
star of ml =0.5 M0, the number of weak transfer events is approximately three times smaller in each case 
(because the weak capture probability in this case is 5-10-4 instead of 1.5-10-3 ~ see Table 1). 
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Note that the results in Table 2 refer to the weak transfer of solid material between neighboring 
planetary systems, but it does not account for the probability of landing on a terrestrial planet in the target 
system. Melosh (2003) and Adams & Spergel (2005) estimate that the latter probability is ~ 10~^. 



6. Summciry and Discussion 

6.1. Summary of the weak transfer mechanism 

We have explored a mechanism that allows the transfer of solid material between two planetary systems 
embedded in a cluster. This mechanism is based on the chaotic dynamics of the restricted three-body model 
of the meteoroid, the star and the most massive planet in the planetary system {Pq, Pi, P2), in which a 
chaotic layer replaces the regular parabolic trajectories of the two body problem of {Pq, Pi). Similarly, 
there is a chaotic layer around the target star (assumed to harbor a planetary system). Weak transfer 
takes place within these chaotic layers because the trajectories have low escape velocities and therefore the 
capture probability in enhanced. We have applied this mechanism to the problem of planetesimal transfer 
between planetary systems in an open star cluster, where the relative stellar velocities are sufficiently low (^ 
1 km/s) to allow weak escape and capture. We found that weak escape and capture within an open cluster 
can enhance drastically the probability of transfer compared to the scenario described in Melosh (2003) and 
corresponding to the exchange of meteoroids between field stars using h}rperbolic trajectories: while Melosh 
(2003) estimated a cross-section of 0.025 AU^ (for capture by a planetary system with a Jupiter-mass planet 
at 5 AU) , and order-of-magnitude estimate for weak capture cross-section would be tt • iJ^^p ~ 10^ AU^ (i.e. 
many order of magnitudes larger) . 

To obtain quantitative estimates, we adopt the average cluster properties inferred for the Solar birth 
cluster (Adams 2010), with N = 4300 members, a total mass of Mduster = 3784 Mq (using an average stellar 
mass of 0.88 Mq resulting from the initial mass function) and a cluster scale length of a = lpc^N/300 = 
3.78 pc. Such a cluster is expected to have a lifetime of 2.3 (Mciuster/A/©)"'^ Myr = 322.5 Myr (ranging 
from 135-535 Myr, for N = 1000 10000). 

With the aid of Monte Carlo simulations, we estimate the probability of weak capture for meteoroids 
that have weakly escaped their original planetary system and are weakly captured at a neighboring system. 
We consider three cases where the source Pi and the target star P* have masses of mi = 1, m* = 1 (Case 
1), TOi — 1, ml = 0.5 (Case 2) and mi = 0.5, ml = 1 (Case 3). The resulting weak capture probabilities arc 
0.15%, 0.05% and 0.12% respectively. This capture probability is much larger than the capture probabilities 
obtained in previous studies; for example, Adams & Spergel (2005) found capture probabilities of 10~^ 
for mean ejection speeds of ~ 5 km/s, typical of hyperbolic ejecta of the solar system. 

Adopting parameters from the minimum mass solar nebula (Weidenschilling 1977; Hayashi 1981), 

considering a range of planetesimal size distributions derived from observations of asteroids and KBOs and 
theoretical coagulation models (Bernstein et al. 2004; Bottke et al. 2005; Kcnyon et al. 2008), and taking into 
account the results from Oort Cloud formation models (Brasser et al. 2011) for the fraction of planetesimals 
that are subject to weak escape from the early solar system, we estimated the number of meteoroids that 
may have been delivered to the weak stability boundary of the solar system over the lifetime of the Sun's 
birth cluster. Using this number and the probability of the weak capture of meteoroids on weakly escaping 
orbits, we calculated the number of weak transfer events from the early solar system to the nearest star in 
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the cluster, assuming it is a solar-type and harbors a planetary system. This number depends strongly on 
the adopted planetesimal size distribution. We find that, for the cases where the power law size distribution 
at the small sizes (dN/dD ex. for D < Dq) has index q2 = 3.5, the expected number of weak transfer 

events between two solar type stars is of the order of 10^'*-3 • 10"'^^; for a shallow size distribution (52 = 1-1) 
the number is of the order of 10^. 

We conclude that solid material could have been transferred in significant quantity from the solar 
system to other solar-type stars in its birth cluster via the weak transfer mechanism described here. 



6.2. Implications for lithopanspermia 

Section |6.1| indicates that the transfer of planetesimals via the weak transfer mechanism described in 
this paper is likely to be the dominant process for the exchange of solid material amongst planetary systems 
in a star cluster. This is of interest for lithopanspermia because if life arose in any of these systems before the 
cluster dispersed, this mechanism may have allowed the exchange of life-bearing rocks amongst the planetary 
systems in the cluster. Within the context of the solar system's formation and dynamical history, in this 
section we discuss how much material originating from the Earth's crust may have been available for weak 



transfer (Section 6.2.1), and, given the time constraints imposed by the weak transfer mechanism, whether 



or not there was a "window of opportunity" for lithopanspermia to take place (Section 6.2.2) 



6.2.1. Earth crustal material available for weak transfer 

We first obtain an order-of-magnitude estimate of how much material may have been ejected from 
the Earth's crust as a consequence of the heavy bombardment that took place before the cluster dispersed. 
Following Adams and Spergel (2005), we assume that I km of the Earth surface was ejected, with a total 
mass of Mb ^ S^^^^^M^ ~ 5 • 10~^ • l(km)M(^ ~ 3 • 10^"' • l{km) g. Adopting a power-law distribution, 
dN/dm cx m^", with a = 5/3, mi — 10 kg and m2 — 10~^ (corresponding to objects 10 km in size), 
this total mass would be distributed in Nn ~ a~i^2-o, ~ 2 • 10^^ • Z(km) bodies. A significant fraction 

of these fragments would have been ejected on hyperbolic orbits, i.e. beyond the domain of weak escape. To 
estimate how many of these bodies may have populated the weak stability boundary, we use the Oort Cloud 
formation efficiency of ^ 1% (Brasser et al. 2011): of the order of 2 • 10^'^ • Z(km) bodies with a terrestrial 
origin that may have been subject to weak escape. 

An additional factor to consider is that a significant fraction of the Earth ejecta resulting from the 
bombardment would have been heated by shocks to pressures and temperatures high enough to sterilize 
the fragments (up to 50 GPa and several 100 °C respectively). However, a few percent of the material that 
originated from the spall region of the impacts -located a few projectile diameters away from the impact point 
in an area where the shocks cancel out- would have remained weakly shocked, achieving a peak temperature 
< 100 °C that would allow microorganisms to survive (Artemieva & Ivanov 2004; Fritz et al. 2005, Pierazzo 
& Chyba 1999, 2006). In fact, laboratory experiments have confirmed that several microorganisms embedded 
in martian-like rocks have survived under shock pressures similar to those suffered by martian meteorites 
upon impact ejection (StofFler et al. 2007; Horneck et al. 2008). Other laboratory experiments confirm that 
bacteria and yeast spores and microorganisms in a liquid can survive impacts with shock pressures of the 



- 18 - 



order of GPas (Burchell et al. 2004; Willis et al. 2006; Hazell et al., 2010; Meyer et al. 2011). Assuming 
only 1% of the ejected Earth material remained weakly shocked, and factoring this into our estimate above, 
we get that ~ 2 • 10"'^^ • ^(km) life-bearing rocks with an Earth origin may have been subject to weak escape. 



Using the weak capture probability, 1.5 • 10^ , derived in Section 4.4 we estimate that the total number 
of lithopanspermia events between the Earth and the nearest solar-type star in the cluster may have been 
of the order of 3 • 10® • Z(km), where / is the depth of the Earth's crust in km that was ejected during the 
"window of opportunity" . 



6.2.2. Time constraints 

Now it is necessary to discuss the time constrains. We focus on two key aspects: (a) whether there is 
evidence that life may have arisen on Earth before the cluster dispersed, and (b) the survival of life to the 
hazards of outer space during the timescales relevant to weak transfer. 



Is it possible that life arose on Earth before the cluster dispersed? 

The age of the solar system can be established from the dating of its oldest solids, the CAI inclusions in 
C-chondrites, that formed 4.570 ± 0.002 Ga when the refractory elements in the solar nebula first started to 
condense at temperatures of approximately 2000 K (Lugmair & Shukolyukov 2001). Let us assume that the 
stellar cluster was also born at that time. From Hf-W chronometry, it is estimated that the crystallization of 
the lunar magma oceans took place 4.527±0.010 Ga (Kleine et al. 2005), setting the time of the giant collision 
of a Mars-sized proto-planet with Earth that stripped part of its mantle and formed the Moon (Canup 2004). 
The detrital zircons found in Jack Hills in Western Australia show evidence that the Earth may have cooled 
down from this Moon-forming collision when the solar system was ^ 70 Myr old; this evidence comes from 
the heterogeneity of the zircons Hf isotope ratio, ^^^Hf/^^''Hf, a tracer of the crust/mantle differentiation 
(Harrison et al. 2005). Furthermore, the high oxygen isotope ratio, ^®0/^^0, of 3.91-4.28 Gyr old Jack 
Hills detrital zircons suggests that the original rocks formed from magma containing recycled continental 
crust that had interacted with water near the surface. This indicates that liquid water was circulating in 
the upper crust of the Earth when the solar system was only 288 Myr old (Mojzsis et al 2001). Another 
study showed high oxygen isotope ratios in 4.404 Gyr old zircons, suggesting liquid water may have been 
present at an even earlier time, when the solar system was 164 Myr old (Wilde et al. 2001). The temperate 
conditions and possible presence of a hydrosphere indicate habitable conditions that may have allowed life 
to emerge during this period. 

The carbon isotopic ratio of tiny inclusions of graphite in 3.85 Gyr old sedimentary rocks in Greenland 
show an increased ^^CO/^"^CO ratio that is indicative of biological activity, implying that life may have 
emerged before the solar system was 718 Myr old (Mojzsis et al 1996). If this age estimate (based on 
uranium-lead dating of zircons) is correct, this means that life may have been extant very shortly after 
the end of the Late Heavy Bombardment (LHB). Rather than abiogenesis taking place during such a short 
period, this favors the hypothesis that life emerged during the Hadean time and survived the LHB. In fact, 
in a study of the degree of thermal metamorphism suffered by the Earth's crust during the LHB, Abramov & 
Mojzsis (2009) concluded that that it is unlikely that the entire crust was fully sterilized and that a microbial 
biosphere, if it existed, likely survived the LHB. 
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In this section on the implications of weak transfer on the possibihty of lithopanspermia, we will work 
with the hypothesis that life emerged dming the Hadean time. This is motivated by the possible presence of 
oceans under temperate conditions and by the timescales for abiogenesis. Note, however, that evidence for 
life as early as 3.85 Ga as mentioned above is controversial (Moorbath 2005 and references therein). There is 
less controversial evidence of a sulfur-based bacterial ecosystem in Western Australian rocks with an age of 
approximately 3.4 Gyr, i.e. when the solar system was approximately 1170 Myr old (Wacey et al. 2011; see 
also Schopf et al. 1993). But this timeframe would be too late for lithopanspermia via weak transfer to take 
place because at this time the solar maternal cluster would have dispersed. The timescales for abiogenesis 
vary from 0.1-1 Myr for hydrothcrmal conditions at the deep sea, to 0.3 3 Myr for warm puddle conditions 
in shallow water, to 1-10 Myr for subaeric conditions in the soil, at least an order of magnitude less than 
the lifetime of the stellar cluster. 

If life arose on Earth shortly after there is evidence of liquid water on its crust, the "window of 
opportunity" for life-bearing rocks to be transferred to another planetary system in the cluster opens by 
the time liquid water was available, at 164-288 Myr, and ends by the cluster dispersal time, Tduster ~ 
135-535 Myr (Adams 2010). Within this timeframe, there was a mechanism that allowed large quantities of 
rocks to be ejected from the Earth: the ejection of material resulting from the impacts at Earth during the 
heavy bombardment of the inner solar system. This bombardment period lasted from the end of the planet 
accretion phase until the end of the LHB 3.8 Ga, i.e. it finished when the solar system was approximately 770 
Myr old (Tera et al. 1974; Mojzsis et al. 2001; Strom et al. 2005). It represents evidence that planetesimals 
were being cleared from the solar system several hundred million years after planet formation (Strom et 
al. 2005; Tsiganis et al. 2005; Chapman et al. 2007). This period of massive bombardment and planetesimal 
clearing encompassed completely the "window of opportunity" for the transfer of life-bearing rocks discussed 
above and therefore provides a viable ejection mechanism that may have led to weak transfer. 

Survival of life to the hazards of outer space during the timescales relevant to weak transfer 

A final consideration for lithopanspermia is the survival of microorganisms to the hazards of radiation during 

their long journey in outer space. Valtoncn et al. (2009) used a computer model to account for the effects 
of Galactic cosmic rays from all elements up to nickel Z = 28, and for the effects of natural radioactivity 
in meteorites characteristic of Earth and Mars. They found the following maximum total survival times 
in interstellar space (that depend on the size of the parent body): 12-15 M}Tr (for sizes of 0.00-0.03 m), 
15-40 Myr (0.03-0.67 m), 40-70 Myr (0.67-1.00 m), 70-200 Myr (1.00-1.67 m), 200-300 Myr (1.67-2.00 
m), 300-400 Myr (2.00-2.33 m) and 400-500 Myr (2.33-2.67 m). To put these lifetimes into context, note 
that under the best case scenario in which life emerged at the time liquid water was available in the upper 
crust (164 Myr or 270 Myr depending on the authors), lithopanspermia would have had a time window of 
up to 260 -370 Myr (assuming the age of the cluster was in the upper end of ^ 535 Myr). But the survival 
timescales above need to be compared, not to this time window, but to the transfer timescales associated 
with the weak transfer mechanism described in this paper. The latter are as follows: (a) Timescale for 
ejection: the numerical simulations carried by Melosh (2003) indicate that the minimum time between the 
ejection of meteorites from the Earth and exit from the solar system is 4 Myr, with a median time of 50 
Myr. Under our scenario, the time for a meteoroid to exit the solar system can be is estimated using 
Barker's equation, that gives the time of fiight along a parabolic trajectory from periapsis with respect 
to the central star to the distance Resc{m,i). For a periapsis of 5 AU (at the location of the perturbing 
planet) and Resc = l-8-lO^ AU (corresponding to a solar- mass central star), we get that the exit timescale 
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is ^ 6 Myr. (2) Timescale for interstellar transfer: a meteoroid moving at the low velocity of 0.1 km/s 
(typical of the weak transfer mechanism) will take about 3-5 Myr to reach a neighboring star (located at a 
distance D ss 10^ AU). (3) Timescale to land on a terrestrial planet: Assuming that the neighboring star 
also harbors a planetary system, the weak capture mechanism would be active; a captured meteoroid would 
typically need to make multiple periapsis passages, i.e., some tens of millions of years, before collision with 
a planet. We see therefore that the timescales for weak transfer compared to the microorganism survival 
timescales estimated by Valtonen et al. (2009) indicate that the survival of microorganisms could be viable 
via meteorites exceeding ^ 1 m in size. 

It is also of interest to study the exchange of prebiotic molecules between planetary systems, as they 
are more robust to the hazards of outer space. Simple amino acids like glycine have been found in several 
carbonaceous meteorites and in Stardust samples returned from comet Wild-2 (Elsila et al. 2009). Iglesias- 
Groth et al. (2011) argue that amino acids like glycine likely formed in the ISM and in chiral excess, and are 
therefore omnipresenlQ This indicates that they may be available throughout the solar system, increasing 
the volume of material that may be subject to weak transfer. Moreover, there is no reason to assume that the 
fundamental hydrocarbon chemistry from which life developed was not present in other planetary systems at 
the time of their formation. Even though glycine has yet to be detected in the interstellar medium, Lattelais 
et al. (2011) points out that its non-detection (in spite of extensive radio surveys) is probably explained 
because neutral glycine is not the most stable isomer and therefore is not dominant. 

The discussion above assesses the possibility that life on Earth could have been transferred to other 
planetary systems when the Sun was still embedded in its stellar birth cluster. But could life on Earth have 
originated beyond the boundaries of our solar system? Our results indicate that, from the point of view of 
dynamical transport efficiency, life-bearing extra-solar planetesimals could have been delivered to the solar 
system via the weak transfer mechanism if life had a sufficiently early start in other planetary systems, before 
the solar maternal cluster dispersed. An early microbial biosphere, if it existed, likely survived the LHB. 
Thus, both possibilities remain open: that life was "seeded" on Earth by extra-solar planetesimals or that 
terrestrial life was transported to other star systems, via dynamical transport of meteorites. 



APPENDIX A: Earth-to-Moon weak transfer: study case for a large focussing factor 



As noted in Sections 



3.4 



andjij adopting for the focussing factor, Gf = 1 + (^)^, a value of G/ = 2 
is a conservative estimate. This factor can be substantially larger if at a given Rcapi'nT'i), "^oo is smaller than 
the approximate value of U — 1 km/s, while the value of Vesc remains the same as U. Dynamically, the 
way to decrease Vao with respect to P^* , as the meteoroid Pq approaches Pi , is for Pq to decrease its relative 
velocity. This process has been shown to exist, for example, by the operational spacecraft Hiten, that was 
transferred from the Earth to the Moon on a trajectory that goes to ballistic capture at a given distance from 
the Moon. Figure 5 shows a trajectory of the spacecraft leaving the Earth at a periapsis distance of 200 km 



Tglesias-Groth et al. (2011) study is based on the radiolysis and radioracemization rate constants derived 
from laboratory experiments in which glycine was exposed to doses similar to those delivered by the decay 
of natural radionuclides in comets and asteroids during 1 Gyr. The authors extrapolate to the solar system 
age and estimate the original concentration of amino acids at the time of solar system formation, concluding 
that "amino acids were formed in the interstellar medium and in chiral excess and then were incorporated 
in comets and asteroids at the epoch of the Solar System formation." 
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and going to a periapsis distance of 500 km from the Moon after 80 days, where the osculating eccentricity 
with respect to the Moon was 0.945. When it arrived at hinar periapsis, its velocity was approximately Wesc, 
while its Vqo went from a value of 1 km/s to km/s (Belbruno 2004), implying G/ = oo. Figure 6 shows 
the Kepler energy of the spacecraft with respect to the Moon along this transfer. At a sufficiently far 
distance from the Moon, where Voo ~ V^Ek, we see that E). ^ implying Voo ^ 0- 

APPENDIX B: Detailed results from the Monte-Carlo simulations 

For each case, the Monte Carlo simulation ran 5 million trajectories, a sufficient number so that the 
randomization of the initial values produces distributions of the parameters that span their respective ranges. 
The resulting distributions are shown in the histograms of Figures 7-21. Figures 7-11 correspond to Case 1, 
Figures 12-16 to Case 2 and Figures 17-21 to Case 3. The histograms are normalized by the total number 
of cases, so that they represent the probability density function of the parameter being measured (i.e. the 
integral under the curve is 1). Histograms labeled as "all cases" include all trajectories, regardless of whether 
or not capture is achieved, while histograms labeled as "capture conditions" include only those cases that 
end in weak capture near Pf . We now describe some of the features of the histograms corresponding to Case 
1. The other two cases are similar. 

Figure 7 shows that weak capture is inhibited for large initial separations between Pq and Pi (i.e. large 
ko/iD- This is because if escape happens at a larger |ro/i|, it implies a larger v^o and larger relative velocity 
with respect to Pf , decreasing the odds of capture. Note that the input distribution for |ro/i| is non-uniform, 
as we are forcing the kinetic energy with respect to Pi to be positive at t = (which favors higher values of 
ko/il, explaining the shape of the "all cases" histogram); if we were to adopt a uniform initial distribution 
for |ro/i|, we would get a decaying exponential or half-Gaussian distribution for the |ro/i| of the capture 
cases. 

Figure 8 shows the initial velocity of Pq with respect to Pi, corresponding to the velocity distribution 
in Figure 2. [Note that the latter figure shows Voo', to compare the two, the escape velocity (~ 0.2 km/s) 
needs to be subtracted from the abscissa of Figure 8]. The peak in the "capture conditions" histogram 
around 1 km/s is due to the input distributions: the relative velocities of the two primaries were set to be 
less than 1 km/s, while the inertial velocities were forced to be less than 2 km/s. Figure 9 shows that these 
various restrictions cause a peak in the distribution of initial inertial velocities of the primaries at ~ 1 km/s. 
Therefore, the most likely velocity that would allow the particle to approach Pf would be -^1 km/s. The 
distribution of Voo in Figure 2 extends out to 5 km/s; these larger velocities correspond to the small peak 
in the "all cases" histogram at 2 km/s (in Figure 8); however, capture conditions become more unlikely for 
these higher velocities and therefore this peak is not present in the "capture conditions" histogram. 
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Table 1. Weak capture probability from Monte Carlo simulations" 



Case 


Mass of source 


Mass of target 


Number of 


Weak Capture 




star mi (M0) 


star ml (Mq) 


trajectories 


Probability 


1 


1 


1 


5-106 


0.15% 


2 


1 


0.5 


5-10^ 


0.05% 


3 


0.5 


1 


5-106 


0.12% 



*^A11 cases assume a cluster with N = 4300 members and cluster properties 
described in Section 
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Table 2. Estimated number of weak transfer events between the Sun and its closest 

cluster neighbor^ 



Case 




?2 


^0 

(km) 






NwTE(mt = l)d 
(N=4300) 


A 


4.3 


3.5 


100 


1.8-1021 


1.8-1019 


2.7-1016 


B 


3.3 


3.5 


2 


2.7-1020 


2.7-101^ 


4.0-1015 


C 


2.7 


3.5 


2 


8.1-101^ 


8.1-1016 


1.2-101^ 


D 


4.3 


1.1 


100 


5.5-10^ 


5.5-106 


8.2-103 



^Considering four different cases (A, B, C, D) for the planetesi- 
mal size distribution. The size distribution is a broken power-law 
given by dN/dD oc D-«i for D > Do and dN/dD oc D-^^ jf D<Do, 
with Dmax = 2000 km (~ Pluto's size) and Dmin = 1 A*in (~ dust 
blow-out size). 

^^D>26cm is the total number of planetesimals with a diameter 
D > 26 cm , equivalent to masses > 10 kg (assuming p — I 
g/cm~3). 

■^Nr is the expected number of planetesimals with a diameter 
D > 26 cm that populated the primordial Oort Cloud, estimated 
to be ~1% of l^D>2ecm- 

'INwteI'^i = 1) is the number of weak transfer events between 
the Sun and a neighboring solar-type cluster star with ml = 1 
Mq and assuming the cluster has N = 4300 members; this is cal- 
culated multiplying Nr by the weak capture probability in Table 
1 corresponding to = 1 Mq, i.e. Nr x1.5-10~3. 
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FIGURE CAPTIONS 




Fig. 1. Schematic representation of the weak transfer process. It consists on a meteoroid weakly 
escaping from a planetary system, and its subsequent weak capture by a neighboring planetary system in 
the stellar cluster. The meteoroid Pq flybys planet P2 and weakly escapes the central star Pi at a distance 
i?osc('Tii). i?csc(™i) is approximately where the the weak stability boundary of Pi is located, caused by the 
gravitational perturbation of the other A'^ — 1 stars in the cluster. The motion in this region is chaotic and 
lies in the transition between capture and escape from Pi. The meteoroid Pq is then weakly captured by 
the neighboring cluster star Pj* at a distance -Rcap("T'i)i moving to periapsis r* with respect to Pj* (motion 
projected onto a plane). The numerical computation pieces together trajectories from flyby of Pj to the 
distance Posc("T'i) from Pi using RP3D and from Posc(™i) to the distance Rcapiml) from Pj* . It has been 
demonstrated that the piecing together of solutions of two different three-body problems at WSB{Pi) and 
WSB{Pi) can be done in a well defined and smooth manner (see Belbruno, 2004; Marsden & Ross, 2006). 
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Fig. 2. Velocity distribution of 670 low-velocity test particles escaping from the solar system. The 
numerical model is the planar circular restricted three-body problem of the Sun, Jupiter and a massless 
particle. 
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Fig. 3. The diagonal solid lines plot the spatial scale of the weak stability boundary for escape, 
-Resc(wi) = 2Gmi/a^, as a function of stellar mass, mi, for two values of the escape velocity, a = 0.1 and 
0.3 km/s. The diagonal dashed line plots the weak stability boundary for capture, i?cap("^i) = 2Gm*/U'^, for 
the average relative velocity of stars in the cluster, U = 1 km/s. The horizontal dashed-dotted lines indicate 
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the spatial scale of star clusters, Rduster, while the horizontal dotted lines indicate the mean interstellar 
distance, D, for clusters consisting of N^lOO, N=1000 and N = 4300 members. 
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mi*(Msun) 

Fig. 4. The diagonal dashed lines plot the probability for weak capture of meteoroids, G f(Rcap/ D)'^ , 
as a function of the target stellar mass {ml), for three values of the number of stars in a star cluster. For 
Gf, we take the conservative value of 2. The solid lines plot the probability for weak capture of meteoroids 
by a neighbor star of stellar mass equal to that of the source star (mi = ml). The symbols indicate the 
stellar mass distribution. 
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Fig. 5. Weak transfer of the spacecraft Hiten from the Earth (located at (0,0)) to the Moon (located 
at (-1,0)) via the hmar weak stability boundary. This is an Earth-Moon fixed rotating coordinate system, 
projected onto the Earth-Moon plane. The time of travel is 80 days, 14 hours. The x-axis and y-axis are in 
units of 4x10^ km. 
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Fig. 6. Kepler energy of the spacecraft Hiten with respect to the Moon (in units of km^s"^) as a 
function of time (in units of days) along the trajectory shown in Fig. 5. 
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Fig. 7. Results from Monte Carlo simulations of 5 million trajectories between a star Pi of mass 
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mi = 1 M0 and star of mass mf = 1 Mq (Case 1). The Figure shows the probabihty density function 
(normalized to 1) of the initial separation between Pq and Pi (ro/i = |ro/i| at t = 0). 
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Fig. 8. Results from Monte Carlo simulations of 5 million trajectories between a star Pi of mass 
mi = 1 M0 and star Pj* of mass mj = 1 Mq (Case 1). The Figure shows the probability density function 
(normalized to 1) of the initial velocity of Pq with respect to Pi {vq/i — |vo/i| at t = 0). 
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Fig. 9. Results from Monte Carlo simulations of 5 million trajectories between a star Pi of mass 
mi = 1 Mq and star P* of mass m* = 1 Mq (Case 1). The Figure shows the probability density function 
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(normalized to 1) of the initial inertial velocities of Pi and P* with respect to the cluster center {vi = |vi 
and vl = |vj| at i = 0). 
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Fig. 10. Case 1: Probability density function of the initial orientation of Pq about Pi (0o and 60 are 
the spherical angles). The integral under the curve is 1. 
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Fig. 11. Case 1: Probability density function of the angular separation between initial velocities of 
Pi, P*. The integral under the curve is 1. 
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Fig. 12. Same as Figure 7 but for Case 2, with mi = 1 Mq and ml = 0.5 Mq. 
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Fig. 13. Same as Figure 8 but for Case 2, with mi = 1 Mq and ml = 0.5 Mq. 
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Fig. 14. Same as Figure 9 but for Case 2, with mi = 1 Mq and ml = 0.5 Mq. 




Fig. 15. Same as Figure 10 but for Case 2, witli mi = 1 Mq and mf = 0.5 Mq. 
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Fig. 16. Same as Figure 11 but for Case 2, with mi = 1 Mq and ml = 0.5 Mq 
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Fig. 17. Same as Figure 7 but for Case 3, with mi = 0.5 Mq and m* = 1 Mq. 
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Fig. 18. Same as Figure 8 but for Case 3, with mi = 0.5 M© and ml = 1 M©. 
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Fig. 19. Same as Figure 9 but for Case 3, with mi = 0.5 M© and m* = 1 M©. 
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Fig. 20. Same as Figure 10 but for Case 3, witli TOi = 0.5 Mq and ml = 1 Mq 



0.16 




[Capture Conditions 
All Cases 




20 40 60 



100 120 140 160 180 

(deg) 



Fig. 21. Same as Figure 11 but for Case 3, witli mi = 0.5 Mq and ml = 1 Mq. 



